Critical points and non-Fermi liquids in the underscreened pseudogap Kondo model 
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Numerical Renormalization Group simulations have shown that the underscreened spin-1 Kondo 
impurity model with power-law bath density of states (DOS) p{(jj) oc \u!\^ possesses various 
intermediate-coupling fixed points, including a stable non-Fermi liquid phase. In this paper we 
discuss the corresponding universal low-energy theories, obtain thermodynamic quantities and crit- 
ical exponents by renormalization group analysis together with suitable e-expansions, and compare 
our results with numerical data. Whereas the particle-hole symmetric critical point can be con- 
trolled at weak coupling using a simple generalization of the spin- 1/2 model, we show that the 
stable non-Fermi liquid fixed point must be accessed near strong coupling via a mapping onto an 
effective ferromagnetic ScS = 1/2 model with singular bath DOS with exponent rcff — —r < 0. 
In addition, we consider the particle-hole asymmetric critical fixed point, for which we propose a 
universal field theory involving the crossing between doublet and triplet levels. 
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I. INTRODUCTION 

Quantum impurity models were introduced in the con- 
text of dilute magnetic moments in metals and are based 
on the notion that important aspects of the physics of 
extended systems can be captured by spatially local pro- 
cesses. Indeed, the reduction of a complicated many- 
body problem to the study of a single quantum degree 
of freedom coupled to a simplified environment can be 
an important step in understanding more complex phys- 
ical situations. Impurity models form the basis for the 
investigation of lattice systems containing local moments 
- these turn out to be extremely rich, displaying various 
phenomena such as heavy-electron quantum coherence, 
glassiness, and quantum phase transitions. In addition, 
simpler yet interesting quantum effects that take place 
at the level of a single magnetic moment can also be suc- 
cessfully described using (oversimplified) single impurity 
models pj. The fermionic Kondo effect, i.e., the screening 
of an extra magnetic moment by conduction electrons, is 
certainly the hallmark of this single-impurity approach, 
but other aspects like local non-Fermi liquid behavior and 
boundary quantum phase transitions 2] have also been 
put forward. 

In general, single-impurity models can serve as an ex- 
cellent testing ground for novel techniques and physical 
paradigms that may be helpful in the more complicated 
case of extended or dense impurity systems. In this re- 
spect, recent work based on the so-called pseudogap (or 
gapless) Kondo and Anderson models, in which conduc- 
tion electrons obey a semi-metallic density of states that 
vanishes as p{ijj) oc |a;|'' (r > 0), has shown interest- 
ing promises. Indeed, this problem features a plethora 
of non-trivial properties, which can be elegantly cap- 
tured by coniroZ/ed fermionic renormalization group (RG) 
techniques 0, Q . Powerful non-perturbative Numerical 
Renormalization Group (NRG) simulations 0, were 
able to confirm the success of the analytical approach. 
In particular, recent work 0, |^ has exposed that all 
fixed points appearing in these pseudogap models can 
be accessed perturbatively after identifying the value of 



the DOS exponent r that corresponds to their associated 
lower-critical or upper-critical "dimension" . 

While the fixed points of the pseudogap spin- 1/2 
Kondo and single-orbital Anderson models are now un- 
derstood in this language, here we are interested in cer- 
tain extensions of these models for which numerical data 
is also available: the case of a spin-1 impurity has been 
investigated in the extensive work of Gonzalez-Buxton 
and Ingersent |^, and has shown interesting yet unex- 
plained properties. We note that the physics of spin-1 
Kondo models can be realized in multilevel quantum dot 
systems, and theoretical predictions for transport quan- 
tities reflecting the singular Fermi liquid behavior of un- 
derscreened Kondo spins in a metallic host have been 
put forward @- 

In this paper, our objective is twofold: first, we will 
generalize the theories for the critical fixed points, which 
are already present in the spin- 1/2 pseudogap Kondo 
model, to the S — I situation, following Refs. 0, 0. 
This includes two types of fixed points, one occuring at 
particule-hole (p-h) symmetry and being controlled at 
small r, and the other present away from p-h symmetry 
and captured near r — 1 using a general mapping onto 
an effective interacting resonant level model. Second, we 
describe the physics of the additional stable fixed point 
which emerges at particle-hole symmetry, due to the in- 
stability of the strong-coupling symmetric fixed point for 
S > 1/2. Analyzing the leading relevant perturbation, 
we will be able to capture this non-Fermi phase by a 
mapping onto an effective ferromagnetic spin S'ofr = 1/2 
model with singular exponent rcff — —r < 0. Techni- 
cally, we will perform one-loop perturbativc RG calcula- 
tions together with suitable e-expansions 9] to determine 
thermodynamic quantities; the results compare favorably 
with the numerical data obtained in Rcf. 5. The present 
work also illustrates the broad applicability of the gen- 
eral field-theoretic tools, developed in Refs. H, H, in a 
more complicated setting. 

The bulk of the paper is organized as follows: we start 
by presenting the spin-1 Kondo model with a given pseu- 
dogap density of states (Sec.^, and summarize previous 
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results from NRG calculations fSec. ini|) . In Sec. IIVI we 
focus on the case of particle-hole symmetry. It shows a 
critical point accessible from weak coupling (which is a 
generalization to S* > 1/2 of the critical point analyzed 
by Withoff and Fradkin ^]). We also present an analysis 
and strong-coupling calculations for the stable non- Fermi 
liquid fixed point found in the numerical simulations of 
Ref. J, Sec.^presents a novel field theory that describes 
the quantum phase transition between the free-moment 
and strong-coupling phases in the presence of particle- 
hole asymmetry. In all cases, direct comparison to NRG 
data is made. Technical details on the diagrammatic ap- 
proach can be found in several appendices. We will finally 
conclude the paper and discuss some open questions re- 
lated to the pseudogap Kondo model. 



II. THE MODEL 

The pseudogap Kondo model (for a review, see Ref . 0) 
originates from the question of how a magnetic moment 
reacts to the coupling to an electronic bath which has 
depleted weight at the Fermi level, and is of interest 
to various condensed-matter systems such as d-wave su- 
perconductors, graphite, and semiconductors. The basic 
physical idea is that the possibility of observing a Kondo 
effect is weakened by the lack of low-energy states, which 
results in various quantum phase transitions controlled 
by the strength of the Kondo interaction, the amount 
of particle-hole symmetry breaking, and obviously the 
shape of the electronic density of states. We will model 
the latter by the simplified form 

p{e)=Y,6{e-eu) = NoW9{D^-e'), (1) 

k 

which has an algebraic dependence in energy character- 
ized by the number r > (the standard Kondo model 
corresponds to r = 0). Here D is the half-bandwidth and 
iVo = (r + l)/(2Z?'-+i). 

The Kondo model is characterized by an antifcrromag- 
netic coupling Jb to this fermionic bath, and possibly a 
potential scattering term E^). It reads 

ka aa' 

+ii;„^ct(0)c,(0) (2) 

a 

in standard notation, and the subscript b stands for 
"bare" physical quantities and operators, in view of the 
later RG calculations. For E'q = the model is p-h sym- 
metric, and tunes the amount of p-h asymmetry. In 
the following we will be concerned with the case of a 
spin-1 operator Sb and want to analyze the various fixed 
points of this model. 
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FIG. 1: (color online) Renormalization group flow diagrams 
for the 5 = 1 pseudogap Kondo model, deduced from NRG 
calculations 0. The horizontal axis is the Kondo coupling, 
the vertical axis denotes the amount of p-h symmetry break- 
ing. Depending on the value of the bath exponent r, three 
distinct situations are found. The red lines represent phase 
boundaries, the red dots critical flxed points. For a detailed 
explanation of the fixed points see text. 



III. SUMMARY OF PREVIOUS NUMERICAL 
RESULTS 

The comprehensive numerical study of Gonzalez- 
Buxton and Ingersent @ has exposed that the model 
|2Jl displays several zero-temperature phases, separated by 
boundary quantum phase transitions. Their results are 
summarized in Fig. ^ as a function of Kondo coupling 
Jb, particle-hole asymmetry Eq, and the value of the 
bath exponent r. 

Let us focus on the p-h symmetric case, Eq — 0, first. 
At Jb = we have a free spin 1, dubbed local-moment 
(LM) fixed point, which is stable w.r.t a small Kondo 
coupling Jb for r > 0. Flow towards large coupling 
is possible only when Jb exceeds a critical value Jc - 
this phase transition defines a symmetric critical (SCR) 
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point. Contrary to the S = 1/2 case where a symmet- 
ric strong-coupHng (SSC) fixed point (corresponding to 
Js = cx)) is always attained when Jb > Jc, SSC is gener- 
ically unstable in the S* = 1 model, which instead shows 
an intermediate stable fixed point (NFL). It is located at 
a finite value J* > Jc and corresponds to a non-Fermi 
liquid (NFL) phase. This flow structure does not persist 
for all values of r, since at a value Tmax — 0.27, the SCR 
and NFL fixed points collapse and disappear altogether. 
For r > r,jiax, LM becomes then generically stable [l0| . 

Now we turn to finite p-h asymmetry, i.e., potential 
scattering Eq ^ 0. Then SCR is stable against finite 
Eq only for < r < r* (where r* ~ 0.245), whereas 
NFL is always unstable. In the case r > r* , SCR be- 
comes unstable as well, so that an asymmetric cousin of 
the critical point, ACR, exists and controls the quantum 
phase transition between LM and an asymmetric strong- 
coupling (ASC) fixed point. ACR is associated to finite 
values and Eqc of both the Kondo coupling and the 
potential scattering term. 

The numerical simulations of Ref. have also estab- 
lished the behavior of several observables, such as the 
impurity contributions to the entropy and to the Curie 
susceptibility, for the various fixed points. Whereas the 
trivial cases (LM, SSC, ASC) are readily understood, 
all intermediate-coupling fixed points (SCR, NFL, ACR) 
correspond to fractional values of both the ground state 
degeneracy and spin, and cannot be understood on the 
basis of an independent electron picture. The values ex- 
tracted from the Numerical Renormalization Group data 
are shown on Table I, and in the remainder of the pa- 
per we aim at an analytical calculation of such singular 
properties, which characterize the complicated nature of 
the various ground states occuring in the model. 



Fixed point (J, Eo) 


Entropy 


Curie constant 


LM (0,0) 


In 3 


2/3 


SCR (Jc,0) 


In 3 + ar^ 


2/3-0.7r 


NFL (J*,0) 


In 2 + 2r In 2 


l/4 + 0.4r 


SSC (cx),0) 


In 2 + 2r In 2 


l/4 + r-/8 


ACR {J'„Eoc) 


1.6-0.7(1 -r) 


0.5 - 0.04(1 - r) 


ASC (cx),oo) 


In 2 


1/4 



TABLE I: Numerical values of the impurity contribution to 
the entropy and to the Curie constant, for the various fixed 
points of the 5* = 1 pseudogap Kondo model for < r < 1, 
as determined by NRG 0. For LM, SSC, ASC the quoted 
values are exact. 



IV. PARTICLE-HOLE SYMMETRIC MODEL 

In this section we analyze the particle-hole symmetric 
spin-1 pseudogap Kondo model lj2Jl, with Eq = 0. We 
start with the weak-coupling RG which allows to capture 
the SCR fixed point for small r. Then we consider an 



expansion around the symmetric strong-coupling fixed 
point - this will lead us to an effective ferromagnetic 
Kondo model with singular density of states which we 
use to describe the physics of the stable NFL fixed point. 



A. Weak-coupling analysis of the SCR fixed point 

The study of the regime with small Kondo coupling 
and particle-hole symmetry was done by Withoff and 
Fradkin [i| for S' = 1/2. Perturbative RG is per- 
formed around Jb — 0, i.e., the LM fixed point. The 
tree level scaling dimension of the Kondo coupling is 
dim[JB] = — and the RG flow in the vicinity of the 
LM flxed point can be described by the one-loop beta 
function: 

m = n-f + o{f) (3) 

where j is the renormalized dimensionless Kondo cou- 
pling, defined in Eq. [7| below. The beta function Q 
obviously yields an unstable fixed point at jc — r, which 
corresponds to SCR and controls the transition between 
LM and NFL. The correlation length exponent is 

i=r + 0(r2), (4) 
v 

and r = can be interpreted as lower-critical "dimen- 
sion" for the phase transition. 

These results can be expected to persist for S — 
1, since the beta function for the Kondo problem is 
known 11] to be independent of the impurity spin value 
S. However, as we are interested in computing physical 
observables at SCR, which do differ from their 5' = 1/2 
counterparts, we will proceed with a complete RG analy- 
sis for S* = 1 as well. We only sketch here the intermedi- 
ate results, while details on the derivation can be found 
in Appendix IXI 



1. RG procedure 

The analysis is based on a representation of the spin 1 
by a triplet of Abrikosov fermions, fmS- 

where the Hilbert space constraint 

E = 1 (6) 

m=l,0,-l 

is enforced using a chemical potential A ^ oo. 

To proceed with the RG scheme, we relate the bare 
coupling constant and field, Jb and fmB, to renormalized 
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quantities, j and fm (which bear no index to shorten the 
notation within renormalized perturbation theory): 

NaJB = ^i-'-Zy'Zjj, (7) 

fmB = ^y^/m (8) 

where /i is a renormalization energy scale. These equa- 
tions define renormahzation factors Zj and Zf which 
will be determined using the standard field-theoretic RG 
scheme |14| . employing dimensional regularization and 
minimal subtraction of poles. 

A one-loop calculation (see Appendix of the 
self-energy and the Kondo vertex yields 

Zf = l + Oif), (9) 

Zj = l + ^-+0{f), (10) 
r 

respectively. The RG beta function is obtained by the 
condition that the bare Kondo coupling is scale invariant, 
i.e. dJs/d/i — 0, which gives from {Tj): 

dfi 1+ jj-\Ti{Zj/Zf) 

Inserting (|^- H10|I into this expression proves that the 
beta function in the S — 1 case is still given by Eq. JSJ, as 
expected. We also note that p-h asymmetry is irrelevant 
at this fixed point, with /3(eo) = re^ (where cq cx NqEq 
is the dimensionless potential scattering term). 

2. Impurity entropy 

The impurity contribution to the entropy is obtained 
as the difference of the total entropy and the entropy of 
the electron bath alone. Its T = value represents a 
measure of the impurity ground state degeneracy. For 
the SCR fixed point it turns out that corrections to the 
LM value are extremely small p^ : 

Sf^f^\ni + 0{r^) . (12) 

This makes the comparison to numerical data difficult, 
and we refrain from extracting the prefactor of the 
term. 



3. Impurity susceptibility 

The impurity susceptibility Ximp is similarly calculated 
by removing the free (i.e. Jb = 0) contribution of the 
bulk fermions to the total spin susceptibility (l7| . 

Xtot(T)= fdr{SUr)SUO)) , (13) 
Jo 

with 

^tot - E ^'flBlnB + E ^^4(X)C.(X) . (14) 
m xcr 



Apart from numerical factors, the calculation is analo- 
gous to the S = 1/2 case |T3|, and we find at lowest 
order in j: 

Because we want to obtain the SCR fixed point contri- 
bution, we can replace j = jc = f and take (fi/T)^^ ~ 1 
as r ^ 0, which gives finally the effective Curie constant: 

in agreement with the numerical value found in Tabic I. 

B. Non- Fermi liquid phase near strong coupling 

1. Around the limit J — oo 

To understand the emergence of a stable non-trivial 
fixed point, we consider in the first place the vicinity of 
the SSC fixed point. Indeed, the renormalization fiow 
sketched from the numerical solution demonstrates that 
the strong-coupling fixed point is unstable, which is quite 
curious at first thought. Indeed, at large J, the spin 
1 is partially screened by the fermions at site 0, leav- 
ing a remnant effective spin Ses — 1/2. Applying the 
Nozieres-Blandin argument [3| we see that an effective 
ferromagnetic coupling J^g = —t^^jj is then generated 
between S^n and the electrons cj.(l) sitting on the site 
near cj.(0) (above, toi is the typical hopping between 
cj.(0) and cj.(l), and we assumed a chain representa- 
tion of the electronic bath). Naively we would then ex- 
pect SSC to be stable. This is not the case, because 
at particle-hole symmetry the Green's functions of the 
above electrons are related through: 

(4(0. «^)c^(0, iuj)) ^ =iuj- ioi(4(l. «t^)c<,(l, i^)) 

(17) 

which shows that the fermions c\{\) have an effective 
density of states Pcff(e) oc lej^"*'', with — —r < 0. 

Thus, the effective Hamiltonian near SSC is a ferro- 
magnetic spin-1/2 Kondo model with a negative DOS 
exponent. Interestingly, this model was studied before in 
Ref.IT^, and was found to display a stable intermediate- 
coupling fixed point. To see this, we note that the beta 
function for Jeff reads 

/3(jcff) = ^eff Jeff - Jeff = -rjeff - Jeff • (18) 

This yields a stable fixed point at j*g = —r, which is 
allowed since the effective coupling is ferromagnetic! In 
our original language, this is the NFL fixed point, which 
controls now a whole phase with non-trivial properties, 
and is a unique feature of the 5' > 1/2 pseudogap Kondo 
model. In particular, the above arguments show that 
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this fixed point is located at j* ~ 1/r. The NFL fixed 
point is unstable with respect to particle-hole asymmetry, 
with P{eo) — — reo, a situation which will be discussed in 
Sec. O Finally, we note that the values for r* and r^ax 
obtained by the Numerical Renormalization Group for 
the original problem |^ and for the effective model 
coincide, in agreement with our identification. We would 
like now to compute the physical quantities at the inter- 
mediate NFL fixed point. 

8. Impurity entropy 

We know that the correction to the entropy at a 
particle-hole symmetric Withoff-Fradkin type of fixed 
point is as small as r^. However, we are now expand- 
ing around a strong-coupling fixed point and not with 
respect to a free spin S = 1. Since the impurity entropy 
is defined relative to the decoupled (J = 0) limit, we have 
to take into account the contribution at SSC {Jb = oo), 
which contains a In 2 term due to the underscreened spin 
5cff = 1/2, and depends on r as well ^ due to the 
singular nature of the resonant level limit (although the 
SSC fixed point is trivial). We obtain finally 

^NfpL = In 2 + 2r In 2 + 0(r3) , (19) 
as anticipated numerically (Table I). 

3. Impurity susceptibility 

We have to follow the same argument as above, and 
sum the contribution of the strong-coupling fixed point, 
Txfmp = 1/4 -f r/8 (see 0) as well as the non-trivial 
correction due to the stable 5cff ~ 1/2 ferromagnetic 
fixed point with exponent rcs = —r, which is 
A(Tximp) — —''off/4 — r/A. The complete result thus 
agrees with the value in Table L 

TX^^ = i + I • (20) 

V. PARTICLE-HOLE ASYMMETRIC MODEL 

The nature of the particle-hole asymmetric fixed point 
of the pseudogap Kondo model, present at Ki ^ and 
r > r* , was understood in the recent work [J, Q . The 
crucial point noticed there is that this quantum phase 
transition separates the local-moment fixed point from 
the asymmetric strong-coupling one, as can be checked 
in Fig. n so that the corresponding field theory can be 
constructed by "mixing" the relevant degrees of freedom 
associated to both fixed points. In the 5=1/2 case, 
one needed to consider a level crossing of a doublet (as- 
sociated to LM) and a singlet (from to ASC), with tran- 
sitions allowed through the coupling to the conduction 
electrons. Not surprisingly, the effective theory took the 



form of an infinite-/7, i.e., maximally p-h asymmetric, 
Anderson model. This was perfectly consistent with nu- 
merical data which indicate the phase transitions in the 
pseudogap Kondo and Anderson models are in the same 
universality class. The theory could be analyzed by ex- 
panding around the level-crossing (valence fluctuation) 
fixed point, and it was shown that r = 1 plays the role 
of an upper-critical "dimension" . 

In the following, we want to generalize the idea of 
Refs.0,|^to the spin-1 case. 

A. Derivation of the critical field theory 

In the present problem, the spin 5 = 1 will be de- 
scribed by a triplet of fermions near LM. On the 
other hand, ASC shows underscreening of the spin, and 
gives a remnant spin-1/2 moment, that we can describe 
by a doublet of bosons b^. Having now the important 
variables to describe the transition at ACR, we need to 
write down the effective theory. Several requirements 
have to be imposed: the Hamiltonian should obey spin 
rotation invariance and should reduce to a spin 1 (resp. 
1/2) Kondo Hamiltonian near the LM (resp. ASC) limit. 
One naturally arrives at the foUowinggeneralization of 
the infinite-C/ Anderson model of Ref. 0: 

n = E^fc4.Cfe. + ^/E/™/™ (21) 

ka rn 

+ Vb [f^oshB Ct (0) + fosb^B Ci (0)] + h.c. 

+ V2Vb [flgb^B Ct(0) + fUshB Ci(0)] + h.c. . 

Here, £/ is the parameter used to tune the system through 
the phase transition; at tree level the transition occurs at 
£/ = 0, and corresponds to the crossing of doublet and 
triplet states, which we call the valence fluctuation (VFl) 
fixed point. The restriction to the physical Hilbert space 
is implemented with the constraint 

Q - E f^BfruB + E f^lB^B = 1 • (22) 

One can check easily via a Schrieffer- Wolff transforma- 
tion that H (|21|l evolves onto a spin 1 (resp. 1/2) Kondo 
model when e/ is large and negative (resp. positive), as 
it should be. 

The effective Hamiltonian (|21|l seems to have nothing 
to do with the original model and in particular it is 
not as simply understandable as in the 5=1/2 case con- 
sidered in 6]. However, we will show that it correctly de- 
scribes the critical properties at ACR. A number of facts 
are encouraging: Power counting at the £/ = Vb = 
fixed point shows that the coupling constant Vb has a 
tree-level scaling dimension dim[VB] = f = (1 — r)/2, i.e., 
it is relevant for r < 1, which indicates the possibility of 
a non-trivial critical point in this regime. In contrast, Vb 
is irrelevant and will flow to zero for r > 1, which estab- 
lishes the role of r = 1 as upper-critical dimension where 
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Vb is marginal. Similar to the situation in the spin-1/2 
model, the transition will become a level crossing with 
perturbative corrections for r > 1. For Vb ^ we have 
strong "valence" fluctuations between the doublet and 
triplet levels, leading to an entropy of In 5, in agreement 
with the numerical result shown in Table I. 

Motivated by this we proceed with an RG expansion 
around the decoupled Vb — fixed point. As demon- 
strated below, we find a flow diagram identical in struc- 
ture to the one of the spin-1/2 model. Fig. 1 of Ref. |^ 
with the difference that the VFl fixed point now corre- 
sponds to a crossing of doublet and triplet levels. This 
flow diagram is very similar to the one of a standard (j>'^ 
theory. In particular, the non-trivial critical point occur- 
ing for r < 1 (which can be understood as the analogue of 
the Wilson-Fisher fixed point) will be identified with the 
asymmetric critical point (ACR) of the original Kondo 
Hamiltonian. 



B. RG calculation for r close to 1 

1. RG procedure 

The renormalization method developed in Ref. 6 dif- 
fers slightly from the technique used in Sec. IIV Al and 
Appendix We still need to introduce renormalization 
factors for the hybridization strength and fields: 



gives: 



NoVb 

fmB 
bcrB 



7I/2 . 
^ f Jm J 

1/2 



1/2 ;7-l/2 



as well as mass renormalization terms. 



(23) 
(24) 
(25) 



(26) 



already written using renormalized quantities. We have 
introduced again a running renormalization scale /i and 
used r = (1 — 7')/2. 

The RG will be performed at criticality, i.e., the 
are chosen such to cancel the real parts of the bare self- 
energies. Including the contribution of the counter-terms 
associated to the above renormalization constants, we 
find (see Appendix the /-electron self-energy contri- 
bution at one-loop: 



D- 



lUJ 

2f 



D 



(27) 

By definition, the field renormalization is determined 
by the condition {d/dLu)Tjf{Q) = 0, and criticality is en- 
sured by the vanishing of the mass, S/(0) = 0, which 





= 1 

r 


Zf 


5\f 





(28) 
(29) 



The bosonic self-energy is similarly calculated (only nu- 
merical prefactors differ) and leads to: 



Zb 



(30) 
(31) 



At this order, the hybridization is left unrenormalized, 
Zv = 1. Finally, imposing the scale invariance condition 
dVe/d/i = 0, we find the beta function: 



dv —TV 

T^' i + v±HZvZ-"^z-"') 

fv+^v"^ + 0{v'^). (32) 



This gives the fixed point value {v*)"^ = (2/5)f. The 
correlation length exponent of the LM-ASC transition 
can be obtained similar to Ref. ^, with the result 



-=r-|-0(H) (r<l) 

V 



(33) 



whereas a level crossing (formally = 1) occurs for r > 1. 
The RG flow is thus identical to Fig. 1 of Ref. B 



2. Impurity entropy 

The correction at order to the free energy is evalu- 
ated in Appendix IbI and reads: 



5 



de — tanh — (34) 
e 2T ^ ' 



Using AS'imp = -d{AF)/dT, and taking the hmit f ^ 0, 
we obtain finally: 



5-/ = l„5-^(l-.) 

which is a agreement with the result in Table I. 



(35) 



3. Impurity susceptibility 

The computation of the Curie constant is done along 
the same line, but is more involved due to the prolifera- 
tion of Feynman graphs, since the total spin operator 

^tot = ''^ fmB fmB 



E 



(-ir 



(36) 
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mixes all three possible fields in the calculation. We sim- 
ply quote the final result, derived in Appendix IBl 



Ximp 



1 

2T 



lOT 



de — -r 

y2 



sinh e/T 
1 + coshe/T 



1 



6 sinh e/T 



(1 + coshe/T)^ 



(37) S, 



We finally insert the fixed-point value v* and take the 
limit f ^ 0, which gives: 



Tx^CR^l_3-21n2 
^Ximp 2 50 '> 



(38) 



Note that the leading Curie constant above is associated 
to the trivial susceptibility of a magnetic system that 
consists of a triplet and a singlet. Because of numeri- 
cal prefactors, e.g. in the partition function, it is not 
simply given by the expression valid for a triplet state 
only, namely S(S-|-l)/3=2/3, but by the different value 
1/2. Finally, we check that equation H^^8|l coincides with 
the results in Table I and vindicates the usefulness of the 
effective theory (|21|l in describing the ACR fixed point. 



imp i_o 
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FIG. 2: (color online) Comparison of the impurity entropy 
Simp between the analytical estimates at one-loop order (con- 
tinuous lines) and the NRG data of Ref. |^ (dots and squares) 
for the three different intermediate-coupling fixed points oc- 
curing in the 5 = 1 pseudogap Kondo model. 



VI. CONCLUSION 

In this paper, we have studied the spin-1 pseudogap 
Kondo model using perturbative RG together with renor- 
malized perturbation theory. We have shown that all 
non-trivial fixed points are accessible in the vicinity of 
their associated upper-critical or lower-critical dimen- 
sion. In particular, we have proposed a novel field theory 
governing the critical behavior in the presence of particle- 
hole asymmetry, and we also have provided a physical 
picture for the stable particle-hole symmetric non-Fermi 
liquid phase (which is not present in the S" = 1/2 case) 
using a strong-coupling mapping to an effective ferro- 
magnetic spin-1/2 model with singular density of states. 
The generic instability of the symmetric strong-coupling 
fixed point SSC and the emergence of the non-Fermi liq- 
uid phase presents a nice example of a situation where 
the purely local strong-coupling picture of Kondo singlet 
formation is clearly inappropriate. We also note that, 
from the arguments given above, our qualitative results 
are expected to be generic to the pseudogap Kondo model 
for all spin values S greater or equal to 1. 

The collected analytical results for thermodynamic 
quantities established in this paper can finally be com- 
pared to the numerical values obtained from the simula- 
tions 5], and are displayed in Fig. |21 and 13 The agree- 
ment is excellent, even at one-loop order, demonstrating 
the usefulness of e-expansions for impurity models with 
power-law bath spectra. 

Using the present RG methods other observables can 
be calculated as well 0, like the critical exponents for 
the local susceptibility or the conduction electron T ma- 
trix. We only mention here the result for the T matrix: 
similar to the spin-1/2 model the critical exponent can be 
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FIG. 3: (color online) As Fig. |21 but for the effective Curie 
constant Tximp. 

determined exactly to all orders in perturbation theory, 
with the result 

Imr(w)cx— (0<r<l) (39) 

which holds at all intermediate coupling fixed points 
(SCR, NFL, ACR). At the critical dimensions, r = 
and r = 1, logarithmic corrections to this result occur. 

The application of our results to experiments requires 
host systems with pseudogap DOS. Candidates are d- 
wave superconductors or materials with a special semi- 
conducting band structure (as planar graphite) - these 
obey a linearly vanishing density of states. Systems with 
exponent < r < 1 are difficult to realize; proposals have 



8 



been made e.g. using disordered mesoscopic conductors 

la- 

On the theoretical side, the single-impurity pseudo- 
gap Kondo problems are now reasonably well under- 
stood. One open issue is a possible analytical descrip- 
tion of a collapse of intermediate-coupling fixed points, 
which occurs both as r — > r*^ (where two critical fixed 
points merge) and at r ^ ?'max~ (where the stable and 
the critical p-h symmetric fixed points meet and disap- 
pear). One idea comes to mind when comparing the en- 
tropies for a spin-S model: the critical fixed point has 
Simp{Jc) = ln(2S' + 1) + 0{r^), whereas for the stable 
NFL fixed point has 5'imp(J*) = ln(25) -f 2r In 2 -f ©(r^). 
The two values cross at fmax ^ 1/(5 In 8), a value that 
becomes small if 5 ^ 1. One may hope that a controlled 
weak-coupling theory describing the physics of Tmax or 
r* exists for large S, although this remains to be demon- 
strated. Finally, we point out that generalizations of the 
pseudogap model to several impurities can show an inter- 
esting interplay of critical points with different nature. 



Acknowledgments 

We thank L. Fritz and M. Kircan for useful discus- 
sions. This research was supported by the DFG Center 
for Functional Nanostructures and the Virtual Quantum 
Phase Transitions Institute in Karlsruhe. 



APPENDIX A: WEAK-COUPLING RG 

Here we provide calculational details for the weak- 
coupHng RG of Sec HvH 



1. Vertex renormalization at one loop 

We first rewrite the bare Kondo Hamiltonian Q in 
terms of renormalized quantities, allowing to separate 
the bare action Sb into a renormalized Kondo part S, 
a constraint term S\, and a series of counter-terms SS: 



Se 



S + Sx + SS, 



(Al) 



S = / dr ^ Cfc^ {dr + ek)Cf^^ + ^ fmdrfm 

k 

+ /d^E r|^Lv.'/,n/™c.(0)c.(0), (A2) 

mm aa 

Sx = Jdrx(^^f„J„,-?j , (A3) 
5S = [drdXY^Lfm (A4) 
+ / dr {Zj -I) rilv.'/m/mc,(0)c.(0) . 

mm'aa' 
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FIG. 4: (color online). Renormalization of the vertex F^p-'j^^ at 
one-loop level. Continuous and dashed lines denote the free 
Green's functions of the /m and Ca fields respectively, and the 



dot represents the vertex counter-term (Zj — \')V\ 



(0) 



With a renormalized dimensionful coupling J = j/iVo 

.(0) 

lOiT 



the tree level vertices satisfy F[n''it = ri*^, it = Fn^^ 



O-UT ^ OlTi ^ 

;^PIj^ — P J/ a/2 and F^Lo-o- = /x"'' Jto(-1)'^/2, with 
all other terms being zero. Indeed, inserting the defi- 
nitions 0-© into the above expression for the action, 
one recovers the bare action of the problem. The free 
propagators are given by: 



Gfoiiuj) = 



1 



ILU - Ck 
1 



(A5) 
(A6) 



where the limit A — > oo should be taken only at the end 
of the calculation |l3j| . The diagrammatic expansion at 
one loop of the vertex Fioxf given in Fig. 0] 

Evaluating the loop at zero temperature, one finds: 



.(1) 
lOiT 



J 



-2r 



j2 I 

71 //^"V 



V2 



(l + rln(|H/M)) 



(A7) 



where we have used j — No J, and the limit r has 
been taken in the last integral. The renormalization con- 



dition for the vertex is f[q''|^-|.(0) = rio\| — M ^J/V^, so 
that the result ()10|l follows. 



2. Impurity susceptibility 

The diagrammatic expansion of the impurity suscepti- 

111 



bility follows from its definition (|13I14(I . and accounting 
for all numerical prefactors (which arise either from the 
(— 1)'^/2 term in (|14(l or from combinatorial reasons), we 
have to calculate the contribution shown in Fig. Note 
that (Q)ximp and not Ximp has a simple diagrammatic 
expansion, since the constraint Q = flisfmB ~ ^ 
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means that we are not expanding around a free-particle 
limit. At lowest order, we have (Q) = 3exp(— /3A). Com- 
puting the Matsubara sums at finite temperature, one 
gets the result (|15|) in the limit where the half-bandwidth 
D goes to infinity. 



APPENDIX B: RG AROUND THE 
VALENCE-FLUCTUATION LIMIT 



This appendix contains details for the RG calculation 
of SecEU 




FIG. 6: (color online). Self-energy at one loop for the /i 
fermion (upper part) and 6| boson (lower part) respectively. 
Wiggly lines represent the bosonic free propagators. Crosses 
and dots indicate respectively mass and field renormalization 
counter-terms. 



1. Pseudoparticle self-energies at one loop 

Again, we will set up notations and express the bare 
action Sb = S + S\ + SS in terms of renormalized quan- 
tities: 

S = IdrY^^kAdr + ekK, (Bl) 

■Jo m 

ff^ - 

+ UtV ^/[hh,cAQ) + hh.c^iQ)] -fee, 
Jo 

+ dTV2V ^l^[hb^c^{0) + J-ib^cAQ)] +C.C., 
Jo 

5a - /dr Af ^ + ^foA-lV (B2) 
5S = [dT5XfJ2f„J„, + 5XtY,Kb. + --- (B3) 

■^0 m a- 

As above, we have introduced a renormalized dimension- 
ful hybridization V = v/Nq. We did not write the vertex 
renormalization factor, since it will not appear at this 
order of the calculation. We also have to introduce the 
bosonic Green's function [13: 



1 



11/ — A 



(B4) 



2x 




4x 




t 



Calculating the self-energy for the /i fermion shown in 
Fig. 1^ we find at zero temperature: 



+ SXf - iLoiZf - 1) 



Jo e-iuj 

(B5) 

which reduces to H27|) when the limit r — > is taken. 
The calculation for the h-\ proceeds along the same line, 
although the prefactors differ slightly (see Fig.|Hl). 



2. Impurity entropy 



To evaluate the impurity entropy, we need to expand 
the physical partition function which is represented 
diagrammatically in Fig. [7| The free energy can then be 
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FIG. 5: (color online). Diagrammatic expansion of (Q)ximp 
at lowest order for the Kondo model at SCR, including the 
complete prefactors. 



FIG. 7: (color online). Expansion of the partition function Z 
up to order and SXf^t, including the combinatorial factors. 
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which is the result quoted in Eq. (|34(l from which the 
impurity entropy at the ACR fixed point foUows. In de- 
riving this result we have used the renormalized values 
for the counter-terms (|29|I - H31I) . which provide the cor- 
rection in the first line of the above expresion. One can 
check that they actually drop out when the entropy is 
calculated, but those corrections are nevertheless crucial 
for regularizing the forthcoming impurity susceptibility. 



3. Impurity susceptibility 



FIG. 8: (color online). Diagrammatic expansion of (Q)ximp 
at lowest order for the effective theory describing ACR, in- 
cluding the complete prefactors. 



obtained through F = — TlnZ, and reads: 



F 



12 
— 1 
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f^'' d6 




Jo 
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de — tanh 
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2T 



(B6) 



Computing the impurity susceptibility is quite cumber- 
some, since many graphs appear at intermediate stages 
of the calculation, due to the fact that both pseudopar- 
ticles now carry a spin index. First, we recall that the 
quantity which admits a simple diagrammatics is in fact 
(Q)Ximp rather than Ximp, where <^Q) is obtained from 
(Q) — &~^^Z and the above result for the free energy, 
Eq. 1B6II . Again, we simply quote the collection of graphs 
that need to be evaluated in Fig. |S1 (with the correspond- 
ing prefactor), and refer the reader to Ref. for more 
details. 
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